To what systems does the Bohigas conjecture apply? 
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We test the applicability of the Bohigas conjecture to systems whose Hamiltonian is not written as 
a closed form analytic expression. A class of such Hamiltonians is created and appear to violate the 
conjecture. Numerical methods are employed to find the spectra of a two-dimensional, classically 
chaotic "billiard" system whose Hamiltonian is in this class. We find that the spectral fluctuations 
are not in agreement with the conjecture. 
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I. INTRODUCTION 



There has been considerable interest in the subject of 
quantum chaos for more than a quarter century. One of 
the critical questions in the field is what — if anything — 
are the definitive signatures in a quantum spectrum that 
reveal the underlying classical dynamics to be chaotic [lj]. 
The so-called Bohigas conjecture^ has played a pivotal 
role in the field. It proposes that the statistical properties 
of spectral fluctuations are given by random matrix the- 
ory (RMT) provided that the associated classical dynam- 
ics are strongly chaotic. In particular, the quantum level 
statistics for time-reversal invariant Hamiltonians which 
correspond to classically chaotic systems are given by a 
gaussian orthogonal ensemble (GOE). The level statis- 
tics for time-reversal non-invariant systems are given by 
a gaussian unitary ensemble (GUE) Q. A critical fea- 
ture of the level statistics with these ensembles is the 
phenomenon of level repulsion. There are strong correla- 
tions between levels which make it far less likely to find 
two nearly degenerate levels than would be the case if the 
energy levels were uncorrelated with fixed average den- 
sity (Poisson statistics). The study of level statistics as 
a probe into the underlying dynamics is quite old, going 
back to seminal work on nuclear spectra by WignerQ, 
MehtaQ and DysonQ. The underlying assumption in 
this early work was that the statistical analysis was jus- 
tified by the complexity of the underlying system. The 
key insight in the seminal paper of Bohigas, Giannoni, 
and Schmit (BGS) was that even simple systems such 
as two-dimensional billiards should have level statistics 
described by RMT provided the underlying dynamics is 
chaotic and the system had no discrete symmetries. 

Since the publication of the Bohigas conjecture numer- 
ous chaotic systems have been studied and found to fol- 
low the RMT predictions. For example, the energy lev- 
els of hydrogen atoms in strong magnetic fields and the 
energy levels of complex nuclei both follow RMT. Since 
the pioneering work of BGS the level spacings of classi- 
cal two-dimensional billiard systems with strong classi- 
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cal chaos have been extensively studied as a laboratory 
to probe the conjecturejl]. This is largely due to the 
relatively tractable nature of these systems. These nu- 
merical studies have been in concord with the Bohigas 
conjecture. 

Given the overwhelming success of the Bohigas conjec- 
ture, it is widely believed to be true. Over the years there 
have been numerous attempts to prove the conjecture or 
aspects of it from first principles. One strategy is based 
on mapping the chaotic system into a zero-dimensional a- 
model mechanics @- Another is based on the special role 
of periodic orbits in the Gutzwiller trace formula 11| . It 
has been argued that given certain technical assumptions 
this approach is adequate to demonstrate that the conjec- 
ture holds in the sense that certain universal correlations 
in the level density which match RMT can be obtained 
Despite these promising formal developments, it 
should be clear at the outset that chaotic classical dy- 
namics by itself does not cause the quantum system to 
obey RMT level statistics. Classically chaotic systems 
with discrete symmetries obey Poisson rather than RMT 
statistics for an obvious reason: the different symmetry 
classes of the quantum system are uncorrelated. To gain 
insight into the question of what does cause RMT statis- 
tics to emerge from classically chaotic systems it is useful 
to understand the class of systems for which the Bohigas 
conjecture applies. 

It is usually believed that the Bohigas conjecture holds 
far more generally than solely for nonrclativistic quantum 
mechanics with Hamiltonians of the standard form. For 
example, the conjecture has been tested experimentally 
in microwave cavities where the classical mechanics is the 
relativistic motion of photons with specular reflection off 
of walls and the quantum levels are simply the normal 
modes of the cavity [l^. This paper addresses a critical 
aspect of this question, namely, whether the conjecture 
applies to the level statistics of all quantum mechanical 
systems whose associated classical dynamics are strongly 
chaotic and which do not have discrete symmetries. In 
some sense it is already known that it does not. Long ago 
it was shown that chaotic billiards on pseudo-spheres (ge- 
ometries with negative curvature) can have spectral prop- 
erties inconsistent with RMT[13j. However, the impor- 
tant question still remains as to the circumstances that 
the Bohigas conjecture should apply for systems with a 
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flat geometry. 

Typically when considering the BGS conjecture the 
quantum systems under study have Hamiltonians of 
rather simple forms which can be written as closed form 
analytic expressions in terms of the quantum position and 
momentum operators. The associated classical Hamilto- 
nian is obtained by starting with the quantum Hamilto- 
nian and replacing the quantum operators for positions 
and momenta with c-number variables. However, in the 
space of possible quantum systems the set of Hamilto- 
nians which can be written in closed form in terms of 
position and momentum operators is an infinitcsimally 
small fraction of the set of all quantum Hamiltonians. 
The question we address is whether the BGS conjecture 
holds once the restriction that the Hamiltonian must be 
written as a closed form analytic expression is dropped. 
We will show that it does not. We will construct an 
explicit class of counterexamples: quantum mechanical 
systems which do not follow RMT level statistics despite 
having classical limits which have chaotic dynamics with- 
out discrete symmetries. 

Before discussing these in any detail, we note at the the 
outset that the class of counterexamples we construct is 
not particularly profound. Indeed, the construction used 
is quite contrived and is in some essentially trivial. It 
basically exploits the fact that mapping from quantum 
to classical systems is not one-to-one. The strategy is to 
find two distinct quantum Hamiltonians where each cor- 
responds to the same chaotic classical systems. One of 
these quantum systems will presumably satisfy the Bo- 
higas conjecture and, by construction, the other will not. 
Despite the contrived nature of this construction, we be- 
lieve that the existence of these rather simple counterex- 
amples may give some insight into the nature of the BGS 
conjecture. 

This paper is organized as follows: in the next section 
we discuss a rather general class of quantum Hamilotni- 
ans which can be uniquely specified but for which there 
is no closed form analytic expression in terms of posi- 
tion and momentum operators. We show that by taking 
a particular limit, one can construct Hamiltonians with 
spectra that violate the Bohigas conjecture — the under- 
lying classical dynamics are chaotic, it has no discrete 
symmetry, and the system is in the semi-classical regime 
but the statistics of level fluctuations are Poissonian. In 
the following section we provide numerical support that 
this construction does indeed produce Poissionian level 
statistics. Finally we make a few concluding remarks 
about the domain of validity of the Bohigas conjecture. 



II. A GENERAL CONSTRUCTION OF 
SYSTEMS WHICH VIOLATE THE BOHIGAS 
CONJECTURE 

As noted in the introduction, we are interested in 
studying Hamiltonians which are, on the one hand, 
uniquely specified but on the other hand cannot be writ- 



ten analytically in closed form in terms of position and 
momentum operators. One obvious choice for this are 
Hamiltonians expressed in terms of integrals which can- 
not be taken analytically. 

Let us start with some quantum Hamiltonian, Hq 
(throughout this paper a hat will indicate a quantum 
operator). Hq is chosen so that it has no discrete sym- 
metries and is expressible in terms of the position and 
momentum operators. Further, let us suppose that the 
classical Hamiltonian associated with Hq gives rise to 
chaotic dynamics in that all trajectories except an in- 
finitesimally small fraction are chaotic. By assumption 
Hq has a discrete spectrum. We denote the n th eigenstate 
as \n) with associated energy eigenvalue E^'- 

H \n)=E^\n) (1) 

Finally, let us suppose that in accord with the Bohigas 
conjecture the statistics of these eigenenergy level fluctu- 
ations are accurately given by RMT. 

Next, suppose that there exists some additional quan- 
tum operator, O, which, for simplicity, we take to be 
dimensionless. Let us assume that O can be expressed in 
closed form in terms of the position and momentum op- 
erators and that [-ff , 0] 7^ 0. Finally let us assume that 
O is a bounded operator: 

< (i/>\6\4>) < 1 (2) 

for all This assumption of boundedness, while not 
required for the construction of a system which violates 
the Bohigas conjecture, simplifies the analysis. The op- 
erator O is given the Schrodinger representation. From 
it we can obtain a Hciscnbcrg operator (D H (t) which cor- 
responds to its time evolution under the dynamics of Hq 
and satisfies the differential equation 

subject to the initial condition O (0) = O. The formal 
solution to this is O h (t) = exp(iHot/h)(D exp(— iHQt/K). 
Consider the following Hamiltonian operator: 

H(e, T ) = H + ^f AtO H {t) (4) 

which depends on two parameters, e and r, in addition 
to whatever parameters are needed to specify Hq and e. 
H(c,t) can be considered as a Hamiltonian in its own 
right. The matrix elements of H{e,r) in the eigenbasis 
of Hq are give by 

(n\H(e,T)\m) = 6 nm 

he sin ((£i 0) - 4?) T /h) (n\6\m) (5) 
+ (E^-E^)r ■ 
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The preceding construction has a natural classical 
analog. The quantum operator Hq is associated with 
i?Q lass (g,p); q and p are vectors of appropriate dimension 
of the classical position and momenta. Similarly O is as- 
sociated with O class (q, p) . The unitary transformations 
in Eq. (|4]) correspond to time evolution under the dy- 
namics of H . Thus the the classical analog is classical 
time evolution under the classical Hq. One can define 
O class (q, p, t) as the time-evolved operator. It satisfies 
the differential equation 



dO class (q,p,t) _ , , 



dt 



(6) 



subject to the boundary condition O clasa (q, p, 0) = 
O cla,ss (q, p) , where the braces indicate a Poisson bracket. 
Thus the classical analog of H(e,r) is 

H^ s {q,p;e,r) = Ht ss (q,P) + ^ f dt cIass (q, p, t) . 

Since Hre,r) is a well-defined quantum Hamiltonian 
and H ma {q,pe,T) its classical analog, it is legitimate to 
ask whether this system satisfies the Bohigas conjecture. 
For generic choices of Hq, O, e and r we have no reason 
to suspect that the conjecture should fail. 



A. The large r limit 

The interesting issue concerns what happens when the 
parameters are not generic. Of particular interest is what 
happens for large values of r. One approach is to consider 
the the formal limit of r — > oo. The study of this limit 
illustrates many of the basic issues and we will consider it 
next. However, as will be discussed subsequently, taking 
the limit raises some subtleties which will need to be 
addressed. 

Looking at the form of Eq. ([5|) it is easy to see 
that in that a critical quantity in the large t limit is 



sm 



E 



(o) 



E 



(0) 



-/»)/(( 



R(°) _ 



t/%), it goes 

to zero for En" 1 ^ Em and to unity for E n °^ = Em\ 
Since the Hq has no symmetries one expects no degen- 
eracies and thus one expects that 



lim (n\H(e,T)\m) = (E n a) +e(n\6\m))6 n 



If one were to define the Hamiltonian 



H{e) = lim H(e, r) 



(8) 



(9) 



it is apparent that [H(e, r), Hq) = since they share a 
common cigcnbasis. By construction the eigenvalues of 
H[e) are given by 

H{e)\n) = E n \n) with E n = E n 0) + e{n\6\n) (10) 



Since the operators Hq and O are unrelated, gener- 
ically one expects the fluctuations in the spectrum of 
eigenvalues of Hq to be uncorrelated with the fluctua- 
tions in (n\(D\n). Let us denote A to be the average level 
spacing of Hq; since the operator O is bounded, A is also 
the average level spacing for H(e) (assuming that the 
energy bins for averaging are much larger than e). Let 
us denote So to be the characteristic scale of the level- 
to- level fluctuation in (n\0\n) in the vicinity of the n th 
level: 

So = 



N ((n + k + l\0\n + k + l)-(n+ k\0\n + k) 
\ ^ 



N 



(11) 



where TV is the number of levels used in estimating this 
scale. There are two characteristic regimes of interest. 
The first is when eSo <C A. In this case the contribu- 
tion of O to the fluctuations in the spectrum of H(e) is 
small compared to the contribution of Hq and thus the 
statistics of level fluctuations are expected to be essen- 
tially given by RMT. However, in the opposite limit of 
(So 3> A one expects that the fluctuations are domi- 
nated by fluctuations in e(n|C|n) (as these are presumed 
to be uncorrelated with the spectrum of Hq); one ex- 
pects that the statistics of level fluctuations of H(e) will 
be approximately Poissonian. If the underlying classical 
mechanics are chaotic, this constitutes a counterexample 
to the Bohigas conjecture. 

Next consider the classical analog of H(e). By analogy 
to Eq. (HI) it is natural to define i? class (e) as 



H c 



3 (e) = lim H clasB {e,r) 



112) 



A simple argument shows that the dynamics of H class (e) 
are chaotic provided that the dynamics of £fg lass are also 
chaotic. iJ class ( e ) is obtained by integrating £) class over 
all times subject to the constraint that at t = the 
trajectory passes through q, p. The basic point is that 
the classical dynamics are ergodic and thus any trajec- 
tory comes arbitrarily close to any point in the accessible 
phase space. Thus, when averaging over long times the 
time-averaged value of class should not depend on the 
choice of value of q,p at t = 0; for any initial condition 
one can always change integration variables so that t = 
corresponds to a point in phase space arbitrarily close to 
the chosen one. Since the time averaged class is inde- 
pendent of q, p it does not contribute to the equations of 
motion; the trajectories are thus the same for JJ class ( e ) 
as for Hq. 

Since the system appears to have both chaotic clas- 
sical dynamics and Poissonian quantum level statistics 
it would seem that this violates the Bohigas conjecture. 
Unfortunately, the construction depends on the r — > oo 
limit; this limit is somewhat subtle and raises questions 
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about the nature of the semi-classical limit. This issue 
will be addressed in the next subsection 



B. The large r limit and the semiclassical regime 

When a system is semiclassical one generically expects 
the dynamics generated by the Heisenberg equations of 
motion to match the dynamics of the classical system 
generated by the classical equations of motions — up to 
small quantum corrections. It is for this reason that the 
classical analog of H(c,t) is given by iJ class (e, r). 

It is worth considering what it means for the quantum 
dynamics to match the classical dynamics. In effect, it 
means that in the semiclassical regime the macroscopic 
motion of generic quantum wave packets — as given by 
the expectation value of the position and momenta — is 
described to high accuracy by the classical equations of 
motion. Of course, by the uncertainty principle there 
must be a spread in position and momentum but this 
spread is assumed to be small compared to the classi- 
cal scales of interest. Similarly, the wave packet must 
have a spread in energy — motion in quantum mechanics 
is driven entirely by the phase evolution of components 
with different energies. Again, it is assumed that this 
energy spread is small compared to the classical scales of 
interest. 

Of course, the macroscopic motion of a wave packet 
is not given exactly by the classical equations — there are 
quantum corrections. These are parametrically small in 
the semiclassical regime; they go as Ti to some power. 
However, over time the cumulative effect of these small 
corrections will eventually become important. One 
knows that this will necessarily occur on time scales long 
enough so that the internal dynamics of the packet be- 
comes relevant; i.e., when the time is long enough to 
resolve the phase evolution of different components of 
the packet. Thus for a quantum system, the notion of 
semiclassical motion is problematic on time scales longer 
than K/A, where A is the typical level spacing for that 
system. 

The construction used above took ar^oo limit and 
thus seems to depend precisely on the long-time regime 
where the semiclassical identification breaks down. This 
raises a potentially important question: should H class (e) 
really be regarded the classical limit of H(e)7 Fortu- 
nately, the answer is yes. The reason is slightly subtle. 

The semiclassical limit is essentially the limit of suf- 
ficiently large quantum numbers. For a non-integrable 
system this effectively means sufficiently high excitation 
energy. The issue is what constitutes sufficiently high. 
A typical criterion is that the excitation energy must 
be much larger than A the typical level spacing. How- 
ever, for the case of H(e) the semiclassical regime is only 
reached when the excitation energy is sufficiently high 
that it is much larger than cSq. The key point is that 
once this regime is achieved, the contributions due to 
fluctuations in the time average of (n\(D(t)\n) are small 



on the scale of the classical physics, and the question of 
whether the large r limit is allowable becomes irrelevant 
to the classical dynamics. Moreover, it is always possible 
to reach the regime where E 3> eSo '■ by construction the 
operator O is bounded and thus so is 8q. 

The condition under which the system is both in the 
semiclassical regime and has approximately Poissionian 
level statistics is 



E > eSo 3> A 



(13) 



One useful way to impose semiclassical constraints is to 
fix the energy and then look at the limit of small h. 
Generically A will scale like some power of Ti. Weyl's 
formula gives 



1 



A(E) 



d^ p d^a 



6(E-H(q,p)) 



(14) 



where / is the number of classical degrees (i.e., the di- 
mension of the vectors p and q). Thus, parametrically, 
A ~ 7r. The scaling of So with h can in principle de- 
pend on both O and H a . For any given system we will 
characterize it to be a power, pq: 5q ~ W° . Condition 
(TiUl) will be satisfied provided that e is taken to scale 
with h according to e ~ h^~ po and be sufficiently large 
numerically to ensure that e5o S> A. 



III. NUMERICAL EVIDENCE 

The preceding argument suggests that there exist 
quantum systems with no discrete symmetries and which 
have well-defined classical limits describing chaotic dy- 
namics and which, nonetheless, do not have RMT level 
statistics. In this section we numerically demonstrate 
that this occurs. 

The system we study will be a two-dimensional bil- 
liard system for Hq. Two parabolic boundaries are used 
which leave the billiard system without any spatial sym- 
metries, see Fig. [T] Such systems arc known to be classi- 
cally chaotic; over the years, the level statistics of these 
systems have been studied and were found to be given, 
with high accuracy, by the RMT results. 

A useful and commonly studied statistical property in 
the context of level fluctuations is the nearest neighbor 
distribution. The energy levels of a system are listed in 
order, E\, E2, E3 . . ., and the successive energy spacings 
are found: Si = (E i+ i — Ei)/A, where A is the average 
energy spacing taken over by a large number of nearby 
levels. p(s), is defined as the probability density for ob- 
taining a given value s. Note p(s) immediately reflects 
basic features of level fluctuations such as tendencies for 
level repulsion. If nearby levels are completely uncorre- 
lated, the spacing distribution is Poissonian: 



p(s) = e 



(15) 



Intcgrable systems generically follow this distribution. If, 
however, a system has the spacing properties predicted 



FIG. 1: Diagram of the billiard system used for H . The two FIG 3. Expe ctation value of the operator 6 plotted against 
parabolic walls remove spatial symmetries and induce chaotic energy (arbitrary units) 
classical dynamics. 




FIG. 2: Nearest neighbor spacing probability distribution, 
p(s), for Ho. s is measured in units of A, the average energy 
spacing. Solid line represents Wigner's GOE distribution for 
2x2 matrices, a good approximation to larger dimensional 
cases. 



by Random Matrix Theory its distribution exhibits level 
repulsion. For the case of time-reversal invariant Hamil- 
tonians the appropriate ensemble is the gaussian orthog- 
onal ensemble (GOE). While no known closed form ex- 
pression is known for this, Wigner's distribution [J], 
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(16) 



which is the exact result for 2x2 matrices, is an excellent 
numerical approximation Q. 

First, let us look at the standard case where one uses 
Hq as the full Hamiltonian. In this case one expects 
RMT to work well. As expected, it does: sec Fig. [2J In 
obtaining this result numerically, we replaced the infi- 
nite potentials of the billiard with numerically large step 
functions, computed matrix elements of the Hamiltonian, 
truncated and then diagonalized relatively large matri- 
ces. We did a standard numerical test to ensure that our 
levels were numerically stable. As expected, this system 
reproduced the RMT level statistics expected from the 
Bohigas conjecture. 



Pushing the numerics for any single system to obtain 
a very large number of accurate levels needed for a high 
statisitcs study can be numerically intensive. As our pur- 
poses are essentially qualitative in any event, we studied 
several different systems with sufficiently different shapes 
as to act as distinct systems. We improved the statistics 
needed to get the histogram in Fig. [2] by combining these 
together. The parameters used for the various systems 
are the curvatures of the two parabolic walls. We follow 
the same strategy when using the construction outlined 
in the previous section. We need to choose a bounded 
operator O in order to proceed. A useful choice is 



O 



Pi 



(17) 



We plot the expectation value of O against the energy 
in Fig. [3] It is apparent that the fluctuations are inde- 
pendent of energy and are of order unity in terms of h 
counting. We chose e to be V ~EA, where E is the average 
energy of the states used in the numerical analysis, and 
computed the spectra of H(e) using the same numerical 
methods used in the computation of the spectrum of Hq. 
Note that with this identification, e ~ Ti and thus eO is 
also of order h. This ensures that in the classical limit, 
the eO vanishes and the classical Hamiltonian simply be- 
comes iJQ lass , which is known to be chaotic. For this 
system, we again produce a histogram of the spacing of 
nearest levels: see Fig. 2J 

It is apparent that the level statistics are not consistent 
with RMT and are very nearly Poissionian. This system 
thus appears to be a viable counterexample to the Bo- 
higas conjecture: its classical limit is chaotic, it has no 
discrete symmetries and yet does not follow RMT level 
statistics. 



IV. CONCLUSION 

The preceding analysis helps to clarify the Bohigas con- 
jecture's domain of validity. The findings of BGS do not 
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FIG. 4: Nearest neighbor spacing probability distribution for 
H(e). Solid line represents a Poissonian probability distribu- 
tion. 



universally hold when the restriction that the Hamilto- 
nians be closed form analytic expressions is relaxed. The 
class of Hamiltonians constructed here may not be of 
much practical utility, but the existence of such systems 
may shed some light on the connection of chaos to quan- 
tum physics. 

One thing is clear: the underlying classical chaos does 
not by itself drive the quantum level spacing statistics 
into the RMT regime. This was known from the out- 
set in that the Bohigas conjecture was known to apply 
only to systems without discrete symmetries — regardless 
of whether the underlying classical dynamics is chaotic. 
The present class of system is simply another class of 
examples for which chaotic classical dynamics does not 
imply RMT level statistics. 

As noted in the introduction, the class of counterex- 
amples to the Bohigas conjecture discussed here are quite 



contrived. The classical dynamics fixes the average quan- 
tum level density according to the Weyl formula. How- 
ever, there is an infinite class of quantum mechanical 
systems which have the same classical limit. These will 
all have the average level spacing as given by the Weyl 
formula. On the other hand, they need to have same 
fluctuations around this average. The counterexamples 
discussed here are designed by fiat to exploit this. 

These counterexamples are in and of themselves con- 
trived and of little practical interest. However, they do 
serve to focus attention on an important theoretical issue: 
namely, the domain of validity of the Bohigas conjecture. 
It seems very likely that the conjecture is valid for some 
class of quantum system: i.e., there is some well-defined 
class of quantum system (with chaotic classical dynam- 
ics) for which RMT correctly describes the statistics of 
quantum level fluctuations. The issue is precisely how 
to specify this class. Clearly this excludes systems with 
discrete symmetries and also systems such as those dis- 
cussed in the previous section. 

It is highly plausible that the class of systems for which 
the Bohigas conjecture holds is extremely large. Indeed, 
it is plausible that it contains all except an infinitesimal 
small fraction of possible quantum Hamiltonians with 
well-defined classical limits (assuming there was a sen- 
sible metric to count this). However, even if this is true, 
there remains the question of how to identify whether a 
given system is or is not in the class. It is particularly 
important to develop a set of criteria which might enable 
one to decide a priori whether a given generic quantum 
system is in this class or not without explicitly calculating 
the spectrum. To the best of our knowledge no general 
set of criteria to answer this exist at present. 

TDC thanks the US Department of Energy for support 
via grant no. DE-FG02-93ER-40762. 



[1] See, for example, F. Haake Quantum Signatures of 
Chaos, 2nd Ed. (Springer- Verlag, Berlin, 2000); H.-J. 
Stockmann, Quantum Chaos: An Introduction, (Cam- 
bridge University Press, Cambridge, England, 1999). 

[2] O. Bohigas, M.J. Giannoni, C. Schmit, Phys. Rev. Lett. 
52 (4), 1-4 (1984). 

[3] T. Guhr, A. Muller-Groeling and H.A. Weidenmuller, 
Phys. Rept. 229, 189-425 (1998). 

[4] E.P. Wigner, Proc. Cambridge Phil. Soc. 47, 790 (1951). 

[5] M.L. Mehta, Nucl. Phys. 18, 395 (1960) 

[6] F.J. Dyson, J. Math. Phys, 3, 140 (1962); J. Math. Phys, 
3, 154 (1962); J. Math. Phys, 166, 140 (1962). 

[7] E. Bittner, H. Markum and R. Pullirsch, 
|arXiv:hep-lat/0110222^ 1 

[8] O. Bohigas, Random Matrix Theories and Chaotic Dy- 
namics, Les Houches Summer School Proceedings Session 
52 , ed. by M.-J. Giannoni, A.Voros, and J. Zinn- Justin, 
(North-Holland, Amsterdam, 1989) 

[9] A.V. Andreev, O. Agam, B.D. Simons, and B.L. Alt- 



shuler, Phys. Rev. Lett. 76, 3947 (1996). 

[10] Sebastian Mller, Stefan Heusler, Petr Braun, Fritz Haake, 
and Alexander Altland Phys. Rev. Lett. 93, 014103 
(2004); Phys. Rev. E 72 046207 (2005); Stefan Heusler, 
Sebastian Miiller, Alexander Altland, Petr Braun and 
Fritz Haake, Phys. Rev. Lett. 98 044103 (2007). 

[11] M. C. Gutzwiller, Chaos in Classical and Quantum Me- 
c/iamcs,(Springer, New York, 1990). 

[12] H.-J. Stckmann and J. Stein, Phys. Rev. Lett. 64, 
2215 (1990); H.-D. Graf, H.L. Harney, H. Lengeler, 
C.H. Lewenkopf, C. Rangacharyulu, A. Richter, 
P. Schardt, and H.A. Weidenmuller, Phys. Rev. Lett. 69, 
1296 - 1299 (1992); P. So, S.M. Anlage, E. Ott, and R.N. 
Oerter, Phys. Rev. Lett. 74, 2662 (1995) 

[13] N.L. Balaza and A. Voros, Phys. Rep. 143, 109 (1986); 
Europhys. Lett. 4, 1089 (1987). 

[14] D. Kaufman, I. Kosztin and K. Schulten, Am. J. Phys. 
67(2), 133(1999). 



